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Abstract 


The open-loop characteristics of a Large-Gap Magnetic Suspen- 
sion System ( LGMSS ) were investigated and numerical results are 
presented. The LGMSS considered provides five-degree-of-freedom 
control. The suspended element is a cylinder that contains a core 
composed of permanent magnet material. The magnetic actuators 
are air core electromagnets mounted in a planar array. Configu- 
rations utilizing five, six, seven, and eight electromagnets were in- 
vestigated and all configurations were found to be controllable from 
coil currents and observable from suspended- element positions. Re- 
sults indicate that increasing the number of coils has an insignificant 
effect on mode shapes and frequencies. 


Introduction 

This paper describes the open-loop characteristics of a Large-Gap Magnetic Suspension 
System (LGMSS). The LGMSS is a conceptual design for a ground-based experiment that can 
be used to investigate the technology issues associated with magnetic suspension at large gaps, 
accurate suspended-element control at large gaps, and accurate position sensing at large gaps 
(ref. 1). This technology is applicable to future efforts that range from magnetic suspension 
of wind-tunnel models to advanced spacecraft experiment isolation and pointing systems. The 
LGMSS considered provides five-degree-of-freedom control. The suspended element is a cylinder 
that contains a core composed of permanent magnet material. The magnetic actuators are air 
core electromagnets mounted in a planar array. In reference 2 an LGMSS configuration using 
five electromagnets was described and an analytical model of the configuration was developed. 
This analytical model was used in reference 3 to investigate two LGMSS control approaches. 
In reference 3 the simplifying assumption was made that the change in field and field gradients 
with respect to suspended-element displacements was negligible. The purpose of this paper is 
to investigate the open-loop characteristics of the LGMSS with the change in field and field 
gradients included to determine if they introduce any unusual stability problems. In addition to 
the minimum five-coil configuration, configurations utilizing six, seven, and eight coils are also 
investigated to determine the feasibility of independently controlling the field component that 
determines the magnitude of the highest open-loop frequencies. Each configuration is described, 
and values of field components generated by the electromagnets at the location of the suspended 
core are presented. The method of calculating field components is discussed in the appendix. 

Symbols 

A system matrix (state-space representation) 

B input matrix (state-space representation) 

B magnetic flux density vector, T 

[OB] matrix of field gradients, T/m 

F total force vector on suspended element, N 

F c magnetic force vector on suspended element, N 

F d disturbance force vector on suspended element, N 

F g gravitational force vector on suspended element, N 
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h 

I 

Ic 

i 

KF 

KT 

K zn 

M 

m c 

T 
T c 
T d 
T E 


V 

v 

Wi 

w 2 

X 

X 

x, y, z 
6 

e 

ft 

V 

Subscripts: 

b 

ij 

(ij)k 

max 

o 

x, y, z 


acceleration due to gravity (1 g % 9.81 m/sec 2 ), m/sec 2 

suspension height (suspended-element centroid to top plane of coils), m 

coil current vector, A 

suspended-element transverse moment of inertia, kg-m 2 
imaginary number 

coefficient matrix of field gradient components 

coefficient matrix of field components 

constant representing magnitude of B zn produced by / max 

magnetization vector, A/m 

suspended-element mass, kg 

total torque vector on suspended element, N-m 

magnetic torque vector’ on suspended element, N-m 

disturbance torque vector on suspended element, N-m 

suspended-element rate to Euler rate transformation matrix for a 
3,2, 1 ( 2 , y , x, respectively) rotation sequence 

inertial coordinate to suspended-element coordinate vector-transformation 
matrix 

velocity vector, m/sec 
permanent magnet core volume, m 3 
weighting matrix (eq. (37)) 
modified weighting matrix (eq. (41)) 
state vector for linearized model 
truncated state vector for linearized model 
coordinates in orthogonal axis system, m 
small increment 

Euler orientation for 3, 2, 1 rotation sequence, rad 
angular velocity vector 
gradient operator 

electromagnet axes 

partial derivative of i component in j-direction 
partial derivative of ij partial derivative in fc-direction 
maximum value 
equilibrium condition 

components along x-, y-, 2 -axes, respectively 
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1-8 coil numbers 

Matrix notations: 


[ ] 

matrix 

[ r 1 

inverse of matrix 

{ } 

column vector 

{ } r 

transpose of column vector 

L J 

row vector 

L J T 

transpose of row vector 


Dots over symbols denote derivatives with respect to time; a bar over a symbol indicates that 
it is referenced to suspended-element coordinates. 

Simplified Analytical Model 

T his section presents a simplified analytical model of an LGMSS that is based on the 
model developed in reference 2. Extensions of the model developed in reference 2 include the 
addition of terms related to the change in field and field gradients with respect to suspended- 
element displacements. The equations are simplified by using small-angle assumptions and 
neglecting second-order terms involving suspended-element motion. The suspended element is a 
cylinder that contains a core composed of permanent magnet material assumed to be uniformly 
magnetized along the long axis. The suspended element is levitated over a planar array of 
electromagnets mounted in a circular configuration. Figure 1 is a schematic representation of a 
five-coil system that shows the coordinate systems and initial alignment. The suspended-element 
coordinate system consists of a set of orthogonal x, y. z body-fixed axes that define the motion of 
the suspended element with respect to inertial space. The suspended-element coordinate system 
is initially aligned with an orthogonal x, y, z system fixed in inertial space. A set of orthogonal 
x b -, y b - 1 z b - axes, also fixed in inertial space, define the location of the electromagnet array with 
respect to the x, y, z system. The x h - and y b - axes are parallel to the x- and y- axes, respectively, 
and the z b - and z-axes are aligned. The centers of the two axis systems are separated by the 
distance h. 

Equations of Motion 

The angular acceleration of the suspended element {St} in suspended-element coordinates 
can be written as (see ref. 2) 

{n} = (I) {x} m 

where I c is the suspended-element moment of inertia about the axes of symmetry (y and z), 
and {T} denotes the total torques on the suspended element (Tj- = 0). Since no torques occur 
about the x-axis, % is assumed to be 0 (see ref. 2) and is not included in equation (1). A bar 
over a variable indicates that it is referenced to suspended-element coordinates and { } denotes 
a vector quantity. The torque {T} can be expanded as 

{T} = {T c } + {T d } (2) 

where {T c } denotes the control torques on the suspended element produced by the electromag- 
nets and {T rf } denotes external disturbance torques. The angular rates of the suspended element 
are obtained by integrating equation (1). The suspended-element Euler rates can be written as 

{<?} = [T E ] {«} 
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where T E is the suspended-element rate to Euler rate transformation matrix for a 3, 2, 1 (that 
is, 2 , y, x ) rotation sequence. By using small-angle and rate assumptions, equation (3) reduces 
to 

M = {«> (4) 

where 

T 

{*} = [0 X by b z J (5) 

The translational acceleration of the suspended element {V} in suspended-element coordinates 
can be written as 

{v} = (l/rn c ) {F} (6) 

where m c is the mass^of the suspended element and { F} denotes the total forces on the suspended 
element. The force { F } can be expanded as 

{F} = {F c } + {F d } + {F fl } (7) 

where {F c ( denotes control forces on the suspended element produced by the electromagnets, 
{F d ) denotes external disturbance forces, and {F y } consists of the force acting on the suspended 
element, due to gravity, transformed into suspended-element coordinates. The suspended- 
element translational_rates are obtained by integrating equation (6). The suspended-element 
translational rates {V} in inertial coordinates are given as 

{V} = [T m ]- 1 {V} ( 8) 

where [Tj^j] is the inertial coordinate to suspended-element coordinate vector-transformation 
matrix. By using small-angle and rate assumptions, equation (8) reduces to 

(V) 2 {V} (9) 

where 

{V} T =L£ y z J (10) 


Magnetic Torques and Forces 

The torque on the suspended element in a given orthogonal coordinate system can be 
approximated as 

{Tc} = v ({M} x {B}) (11) 

where v is the volume of the permanent magnet core, {M} is the magnetization of the permanent 
magnet core, and {B} is the flux density at the centroid of the permanent magnet core produced 
by the electromagnets. (See ref. 2.) In suspended-element coordinates, the torque becomes 

{T c }=u({M}x[T m ]{B}) (12) 

Since, as mentioned earlier, the permanent magnet core is assumed to be uniformly magnetized 
along the x-axis (long axis), {M} becomes 

{M} T =LMx o Oj (13) 

Similarly, the force on the suspended element can be approximated as 

{F c } = v ({M} • V) {B} (14) 
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where V is the gradient operator (ref. 2). Taking the dot product and rearranging terms results 

in / ir \ 

{F r } = u[dB]{M} (15) 

where [<9B] is a matrix of the gradients of B. (See ref. 2.) By using the simplified notation 
fij = dfa/dj , [<9B] can be written as 


BxX 

Bxy 

Bxz 

Byx 

Byy 

Byz 

B Z x 

B Z y 

B zz 


(16) 


In the region of the core, from Maxwell’s equations, V x B = 0 which results in Bjj - B Jt . In 
suspended-element coordinates, the force becomes 


{F c } = v [T m ] [dB] [T m ] _1 {M} 


(17) 


By using small-angle assumptions, the relationship — Bj , , and assuming that 0 r 0 (ref. 2), 


the torque and force equations reduce to 

T dg = vM 1 e{-0yB x -B z ) (18) 

T dS = vMz(-9 z B x + B y ) (19) 

F C j> = vMj ( B X x + 2B z Bxy ~~ 2 OyB X z) (^ 9 ) 

F C y = vMx { — 0 Z B X X + B X y -f OzByy — OyByz) ( 21 ) 

F c -z — v M~x (QyBxx + B X z + QzByz — OyB zz ) ( 22 ) 


The torques and forces are functions of rotations and translations about the nominal operating 
point and coil currents {I\,1 2, ■ • • , Bi)\ thus, 


0 z X y z\ T , [I\ h 7 nJ T ) ( 23 ) 

F c J 


Disturbance Torques and Forces 

The only significant disturbance acting on the suspended element is along the 2 -axis and is 
equal to its weight 

{F g } T =L0 0 —m c g\ ( 24 ) 

where g is the acceleration due to gravity. Other disturbance torques and forces are ignored. 
Transforming to suspended-element coordinates gives 


{F^} — [T™] 

Performing the transformation under the assumptions made for equations 

Fgx = °y ( m c9 ) 

Fgy = 0 

Fqz = ~ m c9 


(25) 

(18)-(22) results in 

(26) 

(27) 

(28) 
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Linearized Equations 

The equations of motion are in the form 


{*} = / 


T/Jc 

„ F fm c 



(29) 


where 

{x} r = [n v % e y e z Vs Vy v z x y z j (30) 

The only nonzero disturbance acting on the suspended element is assumed to be {F^}. (See 
eq. (24).) From the previous equations, 


T/Jc 

F/m c 



(31) 


where 


{x} = L Qy 0 Z X y z\ 

(32) 

{I} T =Ul h ... In\ 

(33) 


The equations can be linearized around the nominal operating point X 0 , I G by performing a 
Taylor series expansion. Neglecting second-order terms and subtracting out X 0 results in 


where 


and 


{*x} = M {£X} + [0]{5I} 



Wi = 


fl//c 

0 

0 

0 

0 

0 

0 

0 

0 

. 0 


0 0 0 0 

1/I C 0 0 0 

0 10 0 
0 0 1 o 

0 0 0 1 / m c 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 


0 

0 

0 

0 

0 

l/m c 

0 

0 

0 

0 


0 0 0 O' 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

1 /m c 0 0 0 

0 10 0 

0 0 10 

0 0 0 1 . 


(34) 


(35) 


(36) 


(37) 
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Expanding A results in 


/ = Wi 


which reduces to 


A = Wi 


-dTy/dQy dTy/dttz 
dTj/dtty dTj/dQj 
Otty/dSly 

. dVj/d0.y 


0 0 T y6y TyQ z 0 

0 o T z0 y T z0 z 0 

1 0 0 0 0 

0 10 0 0 

0 0 FxOy ^x8 z 0 

0 0 FyOy FyQ z 0 

0 0 F ldy F z q z 0 

0 0 0 0 1 

0 0 0 0 0 

.0 0 0 0 0 


dTy/dOy dT v /dz - I 


dVj/dzl 


0 

0 

T- 

± yx 

T- 

1 yy 

T- 
1 yz 

0 

0 

T- 

A zx 

T- 

1 zy 

T- 

± zz 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

F- 

F- 

r xy 

F- 
1 xz 

0 

0 

F- 

r yx 

F- 

r yy 

Fyz 

0 

0 

F- 

1 zx 

F- 

r zy 

F- zz 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

1 

0 

0 

0 


(38) 


(39) 


Finally, by using the expressions for torques and forces developed earlier (eqs. (18)-(22) 
and (26)-(28)), A becomes 



'0 

0 

—B x 

0 

0 

0 

0 

B X z 

-Byz 

-B zz ' 


0 

0 

0 

-B x 

0 

0 

0 

B X y 

Byy 

Byz 


1 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

w 2 

0 

0 


2 B X y 

0 

0 

0 

r> 

D '(x x)x 

JD 

JD (xx)y 

B {xx)z 

0 

0 

Byz 

(Byy B XX ) 

0 

0 

0 

rj 

(xy)x 

B(xy)y 

B {xy)z 


0 

0 

( B xx — B zz ) 

Byz 

0 

0 

0 

B(xz)x 

F(xz)y 

B {xz)z 


0 

0 

0 

0 

1 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

1 

0 

0 

0 

0 


.0 

0 

0 

0 

0 

0 

1 

0 

0 

0 
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where the notation dfij/dk — has been used and 


‘ vM x /I c 000 0 0 0 00 0' 

0 vMx/Ic 00 0 0 0 000 

0 0100 0 0000 

0 0010 0 0000 

w 0 0 0 0 vM x /rn c 0 0 0 0 0 

0 0 0 0 0 vM x /m c 0 0 0 0 

0 0 0 0 0 0 vM x /M c 0 0 0 

0 0000 0 0 100 
0 0000 0 0010 

. o 0000 0 000 1. 

Next, expanding B results in 

VdTy/dh dTy/dI 2 ••• dTy/dlnl 

dTyJdh dTj/dI 2 

B = Wj dQy/dh ■ ■ • ( 42 ) 

- dVj/dh 

Expanding the first term in equation (42) gives 

dTy/dh = vM x [■ -0y (dB x /dI \ ) - {dB z /dI x )] (43) 

Because 6 y is small, the term 6 y {dB x /dI\ ) is neglected. Since the fields and gradients are linear 
functions of coil currents, the components of B z produced by coil n of an n-coil system can be 
written as 

Bzn = K zn {In/ 1 inax) (44) 

where / max is the maximum coil current, K zn is a constant that represents the magnitude of 
B zn produced by / ma x> and I n is the coil current. For the total system, B z can be written as 

Bz = (1/^max) {1} (45) 

where 

L^zJ = [K z i Kz 2 Kzn J (46) 

and 

{I} T =Ul h ... In J (47) 

(See ref. 2.) Since the elements of [K 2 J are constants 

dB z /dl = (l/-^max) L K zJ (48) 
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Similar results are obtained for the other fields and gradients. Terms in B related to the identities 
0y = Sly, 9 Z = x = V%, y = Vy, and z = Vj are 0. Then, B becomes 


B = (l//max)W 2 


-|K 2 J 1 
LK,J 
L°J 
Loj 

[K X x\ 

[Kij/J 

[K X z\ 

LOJ 

Loj 

Loj J 


(49) 


Initial Conditions 

The suspended element is assumed to be initially suspended in equilibrium at a distance h 
above the electromagnet array with the suspended-element coordinates initially aligned with the 
inertial coordinates as shown in figure 1. In equilibrium, F x - Fy = 0 and the only force on the 
suspended element is along the 2 -axis and is equal to the suspended-element weight 

Fj = rn c g ( 5 °) 


From equations (18)- (22), we have 

By — B z = B xx = B X y — 0 


(51) 


and 


Brz = 


rnc9 

vM x 


(52) 


In equilibrium, by using the relationship of equation (52), element (5,3) of the matrix in 
equation (40) reduces to 

ZM--2 B xz = -B xz (53) 

vM x 

From equation (45) , the controlled fields and gradients as a function of I 0 can be written as 


(1/ Anax) 


LKyJ ' 
LK.J 
[K xx \ 

l K xy\ 

[K xz \ . 


{Io} — L By B z B xx B X y B xz \ 


(54) 


For a five-coil system, {Io} can be found by direct inversion of the K matrix in equation (54). 
The generalized inverse (ref. 4) is used to calculate {I 0 } for higher order coil systems. This 
produces a solution where the 2-norm of the current vector {Io} is minimized (i.e., minimum 
SI 2 ). If all electromagnets are the same size and have the same characteristics, the resistive 
power is minimized for room-temperature electromagnets and the stored energy is minimized in 
all cases. Note that if the electromagnets are of varying sizes, scaling of the columns of the K 
matrix may be necessary. 

Alternative solutions can be found if different constraints are applied. The solution giving 
the minimum oo-norm of the current vector corresponds to the maximum individual current 
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value being minimized. With fixed core magnetization, this also gives the current distribution 
required for maximum force and moment capability. However, this solution has been found 
to often exhibit discontinuities with varying suspended-element orientation (ref. 5) and was, 
therefore, not addressed in this paper. Once {I„} is determined, the uncontrolled fields and 
gradients required to complete the A matrix can be calculated. 

One of the objectives of LGMSS development is to allow positioning of the suspended element 
through large angles in yaw (6 Z ) up to 360° (ref. 3). As the suspended element is rotated, the 
equilibrium suspension currents will change. In the appendix of reference 2, the equilibrium 
currents are developed as a function of yaw angle and initial torques and forces on the suspended 
element. Rearranging equation (A12) in reference 2 results in 


O/RnaxJ 


[KF1 


where 


L B z By B xx B X y 


By 


rearranging [KT], equation (48) 

x\ - 

LK 


can 


Bxz \ T = (lA'Mj) < 

be written as 


T 

F 


(l//max) 


sin 0 Z (K X J -I- cos 6 Z |_K y J 


[KF] 


{1^} — L By B z B XX B X y B XZ \T 


where 


[KF] 


- cos 


cos 2 9 Z |K XX J 4- 2 cos 6 Z sin 9 Z [K X y J 4- sin 2 9 Z \K yy \ 

n 1 1 1 ( - ? O n \ I TX” I I n f\ It f j 


COS u z |_1V XX J 4" ^ COS v z Sill u z \_rLxy j -r om v z 

s 0 Z sin 6 Z LK. M j -f- (cos 2 0 Z - sin 2 9 Z ) \ K xy \ + cos 0 Z sin 
qos 6 z LK xz J + sinfl^Kj^J 


(55) 


(56) 


(57) 


(58) 


Results and Discussion 

In order to evaluate the open-loop behavior of the LGMSS, the eigenvalues and corresponding 
eigenvectors of the A matrix of equation (34) were calculated for each coil configuration. The 
eigenvalues represent the frequencies of the open-loop modes of the system, and the eigenvectors 
represent the mode shapes (ref. 6). The system was found to be controllable from coil currents 
and observable from suspended-element positions for each coil configuration. 

Five-Coil Configuration 

The value of B xz required for suspension is found to be 0.096 T/m from using equation (52) 
with the parameters given in the appendix. In principle, with the correct coil configuration 
and current distribution, this can be the only electromagnet field or field gradient present. 
Therefore, substituting this value into equation (40) (with all other field and gradient values 
equal to 0) should allow the fundamental modes of the system to be determined. Making the 
substitution and calculating the eigenvalues and corresponding eigenvectors results in the values 
presented in table 1. From the table it can be seen that two fundamental modes occur and both 
involve motion along the x-axis and rotation about the y- axis. Physically, they result from two 
interactions between the displaced suspended element and the B xz gradient. Mode 1, which is 
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an unstable mode, is related to the Ty X term in equation (39). This term initiates a negative 
rotation (pitch) after an ^-displacement. Mode 2, which is a stable oscillatory mode, is related 
to the F xB> term in equation (39). This term initiates a negative ^-displacement following a 
pitching rotation. The mode shapes are shown in figure 2. For convenience, mode 1 is referred 
to as a tumbling mode and mode 2 is referred to as a rocking mode. 

With the value of D xz given previously, the currents required to provide equilibrium 
suspension for the five-coil configuration were found to be 

{V/max} T = L — 0.775i -0.2398 0.6274 0.6274 -0.2398 J (59) 

by using equation (54). These currents were then used to calculate the remaining fields and 
gradients required to determine the elements of equation (40). The results are presented 
in table 2. Note that some values are 0 through symmetry about the £ 2 -plane, some are 
nominally 0 by choice of electromagnet currents, and others are found to be approximately 
0 for this particular configuration. 

Under the assumptions of reference 3, the dominant term in the open-loop equations was B x . 
Adding the value of B x from table 2 to equation (40) and calculating the eigenvalues and 
corresponding eigenvectors results in the values presented in table 3. The frequencies of modes 1 
and 2 from table 1 have been substantially modified by introducing B x , with the result that 
another mode has been added. Mode 3 is an unstable rotational mode about the 2 -axis and is 
caused by the magnetization vector trying to align itself with the applied field (the compass- 
needle effect). Mode 1 represents a combination of this effect and the tumbling motion due 
to B xz . For convenience, modes 1 and 3 will now be referred to as compass-needle modes. With 
only B x present, the frequencies of both compass-needle modes become 12.862 rad/sec, and 
mode 2 is absent. This matches the results in reference 3 and is consistent with the simplifying 
assumptions made therein. The eigenvalues and corresponding eigenvectors of equation (40), 
from using the remaining values in table 2, are presented in table 4. From table 4 the number of 
modes is seen to increase to 5. Mode 5 is a translational divergence in the y-direction because of 
the Fyy term in equation (39). Mode 4 is a stable oscillatory mode in the 2 -direction because of 
the F Jz term in equation (39). The frequency of mode 1, which is the highest frequency mode, is 
increased by a small amount, the frequency of mode 2 is decreased, and the frequency of mode 3 
is unchanged. The mode shapes are shown in figure 2. 

A review of the results presented in tables 1 , 3, and 4 shows that the highest frequency modes 
are compass-needle inodes caused by the presence of B x . The highest frequency of the open- 
loop modes is critical because of the existence of controller time delays, large electromagnet time 
constants, limited frequency response of power supplies, and practical constraints on position 
sensor sample rate. Adding second-order gradient terms to the system model increases the 
number of open-loop modes. However, the additional modes have much lower frequencies and 
do not modify the high-frequency compass-needle terms significantly. Since the highest frequency 
modes arc caused by the presence of B x and since B x is uncontrolled (see eq. (49)) with the 
five-coil system, a decision was made to investigate the control of B x (and therefore the highest 
open-loop frequencies) by increasing the number of coils. The results are presented in later 
sections. 

The equilibrium suspension current can be determined as a function of yaw orientation by 
using equation (57). The results presented in figure 3 show that the variation of current in each 

coil is sinusoidal. Further, since all the coefficients in are identical because of the axial 
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symmetry of the coil configuration, row 2 of equation (57) implies that ^ In — 0. This is also 

n=l 
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true for all other coil configurations investigated herein. The open-loop modes do not change 
significantly as a function of yaw orientation. 

Six-Coil Configuration 

The currents required to provide equilibrium suspension for the six-coil configuration (by 
using eq. (54) and the value of B xz given previously) were found to be 

{I 0 //max} T = L -0.7973 -0.3947 0.3947 0.7973 0.3947 -0.3947 J (60) 

These currents were then used to calculate the remaining fields and gradients required to 
determine the elements of equation (40). The results are presented in table 5. The eigenvalues 
and corresponding eigenvectors of equation (40), by using the values from table 5, are presented 
in table 6. The results are essentially the same as for the five-coil configuration. 

To investigate the control of B x , as discussed previously, K x and B x were added to 
equation (54). Setting B x to 0 and solving for the currents results in 

{I 0 //max} r = L -239.9959 241.1959 -241.1959 239.9959 -241.1959 241.1959 J (61) 

This solution is obviously unacceptable because the magnitudes of the currents are all very large. 
Also, the coils are driven in a +/ — /+ /— sequence which indicates that the field at the suspended 
element is the sum of opposing and nearly cancelling contributions from all coils. This implies 
that independent control of B x and B xz is not a viable approach with this configuration. 

The variation of equilibrium coil currents as a function of yaw orientation, using the 
generalized inverse to solve equation (57), is shown in figure 4. Again, the modes of the system 
do not change significantly as a function of yaw orientation. 

Seven-Coil Configuration 

The currents required to provide equilibrium suspension for the seven-coil configuration (by 
using eq. (54) and the value of B xz given previously) were found to be 

{Io/Zmax} 7, = L -0.8688 -0.5365 0.1897 0.7812 0.7812 0.1897 -0.5365 J (62) 

These currents were then used to calculate the remaining fields and gradients required to 
determine the elements of equation (40). The results are presented in table 7. The eigenvalues 
and corresponding eigenvectors of equation (40), by using the values from table 7, are presented 
in table 8. The results are essentially the same as for the other configurations. 

Solving for the currents required to set B x to 0 as described previously results in 

{W7max} T = |_ -132.6562 119.2667 -82.8908 29.9521 29.9521 -82.8908 119.2667 j (63) 

This solution is also unacceptable because of the high currents required. Again, the coils are 
driven approximately in a +/ — / + /— sequence. (Coils 4 and 5 behave as a single coil.) 

The variation of equilibrium coil currents as a function of yaw orientation is shown in figure 5. 
The modes of the system do not change significantly as a function of yaw orientation. 

Eight-Coil Configuration 

The currents required to provide equilibrium suspension for the eight-coil configuration were 
found to be 


(VW r = [-0.9753 -0.6843 0 0.6843 0.9753 0.6843 0 -0.6843 J (64) 
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These currents were then used to calculate the remaining fields and gradients required to 
determine the elements of equation (40). The results are presented in table 9. The eigenvalues 
and corresponding eigenvectors of equation (40), by using the values from table 9, are presented 
in table 10. The results are essentially the same as for the other configurations. 

Again, solving for the currents required to set B x to 0 results in 

{Io//m*x} T = l -204.016 144.012 0 -144.012 204.016 -144.012 0 144.012 J (65) 

These currents are similar in magnitude to the currents required by the previous configurations, 
and again the coils are driven approximately in a +/ — / + /- sequence. (Coils 3 and 7 have 
zero current.) 

The variation of equilibrium coil currents as a function of yaw orientation is shown in figure 6. 
The modes of the system do not change significantly as a function of yaw orientation. 

Concluding Remarks 

The open-loop characteristics of a Large-Gap Magnetic Suspension System (LGMSS) have 
been investigated and numerical results are presented. The system was found to be controllable 
from coil currents and observable from suspended-element, positions for all coil configurations 
examined. The basic system investigated provides five-degree-of-freedom control of a suspended 
element which is a cylinder that contains a core composed of permanent magnet material. 
The magnetization vector of the permanent magnet core is horizontal (along the long axis 
of the cylinder). The magnetic actuators are air core electromagnets mounted in a planar 
array. The basic system uses five electromagnets since this is the minimum configuration for 
five-degree-of-freedom control. 

The results indicate that the highest frequency open-loop modes are caused by the presence 
of an x-component of the field ( B x ) which is required to produce equilibrium suspension of 
the suspended element. These modes are called compass- needle modes because they result 
from the magnetization vector trying to align itself with the applied field (the compass-needle 
effect). The highest frequency of the open-loop modes is critical because of practical constraints 
such as position-sensor sample rates. Previous analyses have used the simplifying assumption 
that the change in electromagnet field and field gradients with respect to suspended-element 
displacements is negligible. Adding these effects increased the number of open-loop modes from 
three to five. However, the additional modes have much lower frequencies and do not modify the 
high-frequency compass-needle modes significantly. The highest frequency modes, as discussed 
earlier, are caused by the presence of B x . Since B x is uncontrolled with the minimum five-coil 
configuration, it was thought that independent control of B x (and thus the magnitude of the 
highest open-loop frequencies) would be possible by increasing the number of coils. 

Configurations utilizing six, seven, and eight coils were investigated. The results indicated 
that independently controlling B x requires very high coil currents because the solution for 
controlling B x resulted in the coils being driven in a +/ — / + /- sequence. This indicates that 
the field at the suspended element is the sum of opposing and nearly cancelling contributions 
from all coils. The conclusion was made that independent control of B x with a planar array of 
the type investigated and a permanent magnet core with a horizontal magnetization vector is 
not a viable approach. Increasing the number of coils had an insignificant effect on mode shapes 
and frequencies. 


NASA Langley Research Center 
Hampton, VA 23681-0001 
August 18, 1992 
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Appendix 

Electromagnet Fields and Gradients 

This appendix presents, in the form of tables, the components of fields and gradients 
(including second-order gradients) generated by the electromagnets in configurations of five, 
six, seven, and eight coils mounted in planar arrays. It should be noted that the full set of 
components is not included in the tables since = B Jt and B (lJ > k = B( ik ^ = B^ k ^. The 
coils are arranged in a circular configuration and are constrained to fit inside an 8- by 8- by 4-ft 
volume. (See ref. 1.) An example of a five-coil configuration is shown in figure Al. Values of 
fields and gradients are calculated by using the program VF/GFUN (ref. 7), which is a finite- 
element code using integral equation formulation for the solution of the magnetization of iron 
regions. For the results shown in this report, only the pre- and post- processor OPERA was 
required. OPERA calculates the field from complicated conductor geometries through using 
numerical integration techniques. The parameters of the suspended element and core are the 
same as those in reference 8 and are presented in table Al. 

Five-Coil Configuration 

The electromagnet parameters for the baseline five-coil configuration are the same as those 
in reference 8 and are presented in table A2(a). Values of the fields and first-order gradient 
components are given in table A2(b), and second-order gradient components are given in 
table A2(c). A schematic representation of the five-coil configuration is presented in figure Al. 

Six-Coil Configuration 

An extra coil of the same size can be added to the baseline five-coil configuration without 
violating the volume constraint given previously. The parameters for this configuration are 
given in table A3(a). Values of the fields and first-order gradient components are given in 
table A3(b), and values of second-order gradient components are given in table A3(c). A 
schematic representation of the six-coil configuration is presented in figure A2. 

Seven-Coil Configuration 

To add a seventh coil without violating the volume constraint, the coil diameters must be 
reduced. In an attempt to preserve a natural progression of geometry, this reduction was done 
by holding the proportions of the coils and the maximum ampere turns fixed. This resulted in 
an increase in maximum current density. Parameters of the seven-coil configuration are given 
in table A4(a). Values of the field and first-order gradient components are given in table A4(b), 
and values of the second-order gradient components are given in table A4(c). A schematic 
representation of the seven-coil configuration is shown in figure A3. 

Eight-Coil Configuration 

Adding an eighth coil requires a further reduction in coil diameter, which was accomplished 
in the same manner as that used with the seven-coil configuration. Parameters for the eight-coil 
configuration are given in table A5(a). Values of the field and first-order gradient components 
are given in table A5(b), and values of the second-order gradient components are given in 
table A5(c). A schematic representation of the eight-coil configuration is presented in figure A4. 
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Table Al. Suspended-Element and Permanent Magnet Core Parameters 

Core diameter, 0.1016 

Core length, 0.3048 

Suspended-element mass, a m ( >, kg 23.11 

Suspended-element inertia, a I r , kg-m 2 0.6 

Core volume, u, m 3 2.471 x 10 3 

Core magnetization, M Xj A/m 9.5493 x 10 0 

Suspended-element suspension height, h, m 0.9144 

a Mass and inertia quoted include nonmagnetic components surrounding 
the core (ref. 8). 
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Table A2. Five-Coil Configuration 
(a) Electromagnet parameters 


Electromagnet outer radius, rn 0.386 

Electromagnet inner radius, m 0.173 

Electromagnet height, m 0.493 

Location radius, a m 0.7 

Maximum current density, A/cm 2 1535.87 


a Distance from center of array to axis of given coil. 


(b) Electromagnet fields and first-order gradients 


Electromagnet 

B x , 

T 

By, 

T 

T 

B XX, 

T/m 

Bx y, 

T/m 

B X z, 

T/m 

Byy, 

T/m 

« a 

cq EH 

B zz , 

T/m 

i 

0.0216 

0 

-0.0198 

0.0092 

0 

-0.0497 

-0.0306 

0 

0.0215 

2 

.0067 

.0206 



-.0269 

.0118 

.0152 

.0054 

-.0472 



3 

-.0175 

.0127 



-.0046 

-.0191 

.0400 

-.0170 

-.0290 



4 

-.0175 

-.0127 



-.0046 

.0191 

.0400 

-.0170 

.0290 



5 

.0067 

-.0206 



-.0269 

-.0118 

-.0152 

.0054 

.0472 

- 

- 


(c) Electromagnet second-order gradients 


Electromagnet 

^(xx)x^ 

T/m/m 

D {xy)x- 

T/m/m 

B(xz)x' 

T/m/m 

D 

(xy)y' 

T /m /m 

r> 

(xy)z' 

T/m/m 

B (x z)z- 

T/m/m 

i 

-0.0695 

0 

-0.0716 

-0.0573 

0 

0.1276 

2 

-.0498 

-.0452 

.0570 

.011 

-.0416 

.0394 

3 

.0848 

.0058 

-.0222 

.0178 

.0679 

-.1032 

4 

.0848 

-.0058 

-.0222 

.0178 

-.0679 

-.1032 

5 

-.0498 

.0452 

.0570 

.011 

.0416 

.0394 
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Table A3. Six-Coil Configuration 
(a) Electromagnet parameters 


Electromagnet outer radius, m 0.386 

Electromagnet inner radius, m 0.173 

Electromagnet height, rn 0.493 

Location radius, in 0.822 

Maximum current density, A/cm 2 1535.87 


(b) Electromagnet fields and first-order gradients 


Electromagnet 

i B x , 

T 

By, 

T 

B z , 

T 

Bxxi 

T /rn 

B xy, 

T/m 

Bxz, 

T/m 

By y 
T/m 

By Z , 

T/m 

B zz , 

T/m 

i 

0.0201 

0 

-0.0143 

0.0156 

0 

-0.0404 

-0.0242 

0 

0.0087 

2 

.0100 

.0174 



-.0143 

.0173 

-.0200 

.0055 

-.0349 



3 

-.0100 

.0174 



-.0143 

-.0173 

.0200 

.0055 

-.0349 



4 

-.0201 

0 



.0156 

0 

.0404 

-.0242 

0 



5 

-.0100 

-.0174 



-.0143 

.0173 

.0200 

.0055 

.0349 



6 

.0100 

-.0174 

- 


-.0143 

-.0173 

-.0200 

.0055 

.0349 

- 

- 


(c) Electromagnet second-order gradients 


Electromagnet 

O 

T/m/m 

D 

JD (xy)x'> 

T/m/m 

B{xz)x' 

T/m/m 

ID 

(xy): )p 

T/m/m 

B{xy)f 

T/m/m 

B (xz)z> 

T/m/m 

i 

-0.0356 

0 

-0.0775 

-0.0485 

0 

0.0839 

2 

-.0589 

-.0183 

.0175 

.0171 

-.0545 

.0420 

3 

.0589 

-.0183 

.0175 

-.0171 

.0545 

-.0420 

4 

.0356 

0 

-.0775 

.0485 

0 

-.0839 

5 

.0589 

.0183 

.0175 

-.0171 

-.0545 

-.0420 

6 

-.0589 
1 

.0183 

.0175 

.0171 

.0545 

.0420 
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Table A4. Seven-Coil Configuration 
(a) Electromagnet parameters 


Electromagnet outer radius, m 0.358 

Electromagnet inner radius, m 0.160 

Electromagnet height, m 0.457 

Location radius, m 0.883 

Maximum current density, A/cm 2 1785.5 


(b) Electromagnet fields and first-order gradients 


Electromagnet 

B x , 

T 

By, 

T 

B g , 

T 

Bj-.r ■ 

T/m 

Bxy, 

T/m 

BxZi 

T/m 

Byy, 

T/m 

ByZ'f 

T/m 

Bzz, 

T/m 

i 

0.0169 

0 

-0.0103 

0.0162 

0 

-0.0318 

-0.0190 

0 

0.0029 

2 

.0106 

.0132 



-.0054 

.0172 

-.0196 

:0024 

-.0247 



3 

-.0038 

.0165 



-.0173 

-.0076 

.0070 

.0144 

-.0310 



4 

-.0153 

.0073 



.0095 

-.0139 

.0286 

-.0125 

-.0136 



5 

-.0153 

-.0073 



.0095 

.0139 

.0286 

-.0125 

.0136 



6 

-.0038 

-.0165 



-.0173 

.0076 

.0070 

.0144 

.0310 



7 

.0106 

-.0132 


- 

-.0054 

-.0172 

-.0196 

.0024 

.0247 




(c) Electromagnet second-order gradients 


Electromagnet 

^(xx)x' 

T/ m / m 

B(xy)x> 

T/m/m 

B{xz)x' 

T/m/m 

(xy)y' 

T/m/m 

B(xy)z> 

T/m/m 

^(xz)z' 

T/m/m 

i 

-0.0176 

0 

-0.0694 

-0.0400 

0 

0.0570 

2 

-.0499 

-.0002 

-.0048 

.0142 

-.0512 

.0356 

3 

.0254 

-.0340 

.0305 

-.0129 

.0226 

-.0127 

4 

.0333 

.0189 

-.0494 

.0187 

.0414 

-.0514 

5 

.0333 

-.0189 

-.0494 

.0187 

-.0414 

-.0514 

6 

.0254 

.0340 

.0305 

-.0129 

-.0226 

-.0127 

7 

-.0499 

.0002 

-.0048 

.0142 

.0512 

.0356 
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Table A5. Eight-Coil Configuration 
(a) Electromagnet parameters 


Electromagnet outer radius, m 
Electromagnet inner radius, m . 
Electromagnet height, m . . . . 

Location radius, m 

Maximum current density, A/cm 2 


0.339 

0.152 

0.433 

0.952 

1991.83 


(b) Electromagnet fields and first-order gradients 


Electromagnet 

B x , 

T 

By, 

T 

B~, 

T 

Bx: r, 

T/m 

Bxy j 

T/m 

B,Z, 

T/m 

Byy 

T/m 

Byy 

T/m 

Bzz, 

T/m 

i 

0.0144 

0 

-0.0074 

0.0160 

0 

-0.0248 

-0.0150 

0 

-0.0010 

2 

.0102 

.0102 



-.0003 

.0156 

-.0174 

.0003 

-.0174 



3 

0 

.0144 



-.0150 

0 

0 

.0160 

-.0248 



4 

-.0102 

.0102 



-.0003 

-.0156 

.0174 

.0003 

-.0174 



5 

-.0144 

0 



.0160 

0 

.0248 

-.0150 

0 



6 

-.0102 

-.0102 



.0003 

.0156 

.0174 

.0003 

.0174 



7 

0 

-.0144 



-.0150 

0 

0 

.0160 

.0248 



8 

.0102 

-.0102 

> 


.0003 

-.0156 

' -.0174 

.0003 

.0174 




(c) Electromagnet second-order gradients 


Electromagnet 

^(xx)x'! 

T/m/m 

B( X y)x-. 

T/m/m 

B(xz)xi 

T/m/m 

B(xy)y’ 

T/m/m 

T/m/rn 

B(xz)zi 

T/m/m 

1 

-0.0046 

0 

-0.0599 

-0.0326 

0 

0.0367 

2 

-.0363 

.0099 

-.0167 

.0099 

-.0429 

.0259 

3 

0 

-.0326 

.0257 

0 

0 

0 

4 

.0363 

.0099 

-.0167 

-.0099 

.0429 

-.0259 

5 

.0046 

0 

-.0599 

.0326 

0 

-.0367 

6 

.0363 

-.0099 

-.0167 

-.0099 

-.0429 

-.0259 

7 

0 

.0326 

.0257 

0 

0 

0 

8 

-.0363 

-.0099 

-.0167 

.0099 

.0429 

.0259 
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Figure A2. Six-coil configuration. 



Figure A3. Seven-coil configuration. 



Figure A4. Eight-coil configuration. 
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Table 1. Nonzero Eigenvalues and Eigenvectors for Five-Coil Configuration 

With Only B xz Present 



Mode 1 eigenvalues 

Mode 2 eigenvalues 


+/- 7.7995 rad/sec 

+/— 7.7995z rad/sec 

States 

Mode 1 eigenvectors 

Mode 2 eigenvectors 

Qy 

+/+ 1.000 

+/+ 1.000 

ft; 

0 

0 

Oy 

+/- 0.128 

-/+0.128t 

e z 

0 

0 

Vx 

-/- 0.161 

+/+ 0.161 

Vy 

0 

0 

V- z 

0 

0 

X 

-/+ 0.021 

-/+ 0.02H 

y 

0 

0 

z 

0 

0 


Table 2. Uncontrolled Fields and Gradients for Five-Coil Configuration 


B x , T ... . 

Byy, T/m . . 

Byzi T/m 
B zz , T/m . . 
B(xx)x > T/m/m 
B (xx)i r T/m/m 
B(xx)z> T/m/m 
S (xy)x> T/m/m 
£(xy)j,, T/m/m 
T/m/m 
#(x;)x- T/m/m 
B {xz)y’ T/m/m 
#(xz)*, T/m/m 


-0.0419 
- 0.0002 
. . . 0 
. . . 0 
. 0.1842 
. . . 0 
. 0.0003 
. . . 0 
. 0.0615 
. . . 0 
. 0.0003 
. . . 0 
-0.2473 
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Table 3. Nonzero Eigenvalues and Eigenvectors for Five-Coil Configuration 
With Only B xz and B x Present 


States 

Mode 1 eigenvalues 

Mode 2 eigenvalues 

Mode 3 eigenvalues 

+/— 13.596 rad/sec 

+/— \A7Ai rad/sec 

+/— 12.839 rad/sec 

Mode 1 eigenvectors 

Mode 2 eigenvectors 

Mode 3 eigenvectors 

n y 

+/+ 1.000 

+/+ 1.000 

0 

n- z 

0 

0 

+/+ 1.000 

9y 

+/- 0.074 

— /+ 0.223* 

0 

0Z 

0 

0 

+/- 0.078 

Vr 

-/- 0.053 

+/+ 0.490 

0 

Vy 

0 

0 

0 

Vz 

0 

0 

0 

X 

-/+ 0.004 

— / + 0.109* 

0 

y 

0 

0 

0 

z 

0 

0 

0 


Table 4. Nonzero Eigenvalues and Eigenvectors for Five-Coil Configuration 
With All Nonzero Fields and Gradients Present 



Mode 1 

Mode 2 

Mode 3 

Mode 4 

Mode 5 


eigenvalues 

eigenvalues 

eigenvalues 

eigenvalues 

eigenvalues 


+/- 13.669 

+/- 1.794* 

+/- 12.839 

+/— 5.05z 

+/- 2.505 


rad/sec 

rad/sec 

rad /sec 

rad /sec 

rad/sec 


Mode 1 

Mode 2 

Mode 3 

Mode 4 

Mode 5 

States 

eigenvectors 

eigenvectors 

eigenvectors 

eigenvectors 

eigenvectors 

Qy 

+/+ 1.000 

+/+ 1.000 

0 

-/- 0.003 

0 

n- z 

0 

0 

+/+ 1.000 

0 

+/+ 0.0.75 

6y 

+/- 0.073 

-/+ 0.7557i 

0 

+/- 0.001* 

0 

0 Z 

0 

0 

+/- 0.078 

0 

+/- 0.002 

Vi 

— /— 0..758 

+/+ 0.4.75 

0 

1 

1 

p 

o 

o 

0 

Vy 

0 

0 

0 

0 

+/+ 1.000 

Vz 

0 

+/+ 0.001 

0 

+/+ 1.000 

0 

X 

-/+ 0.004 

— / + 0.248z 

0 

0 

0 

y 

0 

0 

0 

0 

+/- 0.399 

2 

0 

0 

0 

-/+ 0.199* 

0 
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Table 5. Uncontrolled Fields and Gradients for Six-Coil Configuration 


B x , T ... . 
Byy , T/m 
B yz , T/m . . 
B zz , T/m . . 
B ( XX)X , T/m/m 
£( xx )y< T/m/m 
%*)*. T/m/m 
T/m/m 

s (xy) 3 /> T/m/m 
£(*»)*> T/m/m 

B(xz)x> T/m/m 
T/m/m 

B {xz)z’ T/m/m 


0.0478 
. . 0 
. . 0 
. . 0 
0.1498 
. . 0 
. . 0 
. . 0 
0.0503 
. . 0 
. . 0 
. . 0 
- 0.2001 


Table 6. Nonzero Eigenvalues and Eigenvectors for Six-Coil Configuration 
With All Nonzero Fields and Gradients Present 



Mode 1 
eigenvalues 

Mode 2 
eigenvalues 

Mode 3 
eigenvalues 

Mode 4 
eigenvalues 

Mode 5 
eigenvalues 


+/- 14.402 
rad/sec 

+/- 1.993i 

rad /sec 

+/- 13.717 
rad / sec 

+/- 4.520?: 
rad/sec 

+/- 2.267 
rad /sec 

States 

Mode 1 
eigenvectors 

Mode 2 
eigenvectors 

Mode 3 
eigenvectors 

Mode 4 
eigenvectors 

Mode 5 
eigenvectors 

Cty 

+/+ 1.000 

+/+ 1.000 

0 

0 

0 

n- z 

0 

0 

+/+ 1.000 

0 

0 

9y 

+/- 0.069 

— /+ 0.502i 

0 

0 

0 

0z 

0 

0 

+/- 0.073 

0 

0 

Vx 

-/- 0.051 

+/+ 0.509 

0 

0 

0 

Vy 

0 

0 

0 

0 

+/+ 1.000 

V- z 

0 

0 

0 

+/+ 1.000 

0 

X 

-/+ 0.004 

-/+ 0.255i 

0 

0 

0 

y 

0 

0 

0 

0 

+/- 0.441 

z 

0 

0 

0 

— /+ 0.221i 

0 
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Table 7. Uncontrolled Fields and Gradients for Seven-Coil Configuration 


B x, T -0.0514 

Byy, T/rn -0.0001 

B yz , T/m 0 

B zz , T/m 0 

B {xx)x< T/m/m 0.1305 

B (xx)y . T/m/m 0 

B (xx)z > T/m/m -0.0002 

5 (zy)x> T/m/m 0 

T/rn/xn 0.0438 

B (xy)z > T/m/m 0 

T/m/m -0.0002 

T/m/m 0 

B {xz)z > T/m/m -0.1728 


Table 8. Nonzero Eigenvalues and Eigenvectors for Seven-Coil Configuration 
With All Nonzero Fields and Gradients Present 



Mode 1 
eigenvalues 

Mode 2 
eigenvalues 

Mode 3 
eigenvalues 

Mode 4 
eigenvalues 

Mode 5 
eigenvalues 


+/- 14.834 
rad/sec 

+/- 2.139* 
rad/sec 

+/- 14.218 
rad/sec 

+/- 4.201* 
rad/sec 

+/- 2.116 
rad /sec 

States 

Mode 1 
eigenvectors 

Mode 2 
eigenvectors 

Mode 3 
eigenvectors 

Mode 4 
eigenvectors 

Mode 5 
eigenvectors 

Vty 

+/+ 1.000 

+/ + 1.000 

0 

+/+ 0.002 

0 

n- z 

0 

0 

+/+ 1.000 

0 

+/+ 0.003 

0y 

+/- 0.067 

-/+ 0.468i 

0 

-/+ 0.001 

0 

e z 

0 

0 

+/- 0.070 

0 

rH 

o 

o 

b 

1 

Vx 

-/- 0.047 

+/+ 0.548 

0 

+/+ 0.001 

0 

Vy 

0 

0 

0 

0 

+/+ 1.000 

Vz 

0 

1 

1 

o 

b 

o 

t— 1 

0 

+/+ 1.000 

0 

X 

-/+ 0.003 

-/+ 0.256* 

0 

0.001 

0 

y 

0 

0 

0 

0 

+/- 0.473 

z 

0 

0 

0 

-/+ 0.238* 

0 
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Table 9. Uncontrolled Fields and Gradients for Eight-Coil Configuration 


B x , T ... . 

Byy, T/m . . 

Byz , T/m . . 
B zz , T/m . . 
B( xx )x, T/m/m 
B (xx )y, T/m/m 
B(xx)z> T/m/m 
B ( X y)x’ T/m/m 
£(xy)ir T/m/m 
£(sy)*. T/m/m 
£(w)x* T/ m /m 
£(**)»> T/m/m 
B (XZ ) Z , T/m/m 


-0.0560 
. . . 0 
. . . 0 
. . . 0 
. 0.1083 
. . . 0 
. . . 0 
. . . 0 
. 0.0365 
. . . 0 
. . . 0 
. . . 0 
-0.1425 


Table 10. Nonzero Eigenvalues and Eigenvectors for Eight-Coil Configuration 
With All Nonzero Fields and Gradients Present 



Mode 1 
eigenvalues 

Mode 2 
eigenvalues 

Mode 3 
eigenvalues 

Mode 4 
eigenvalues 

Mode 5 
eigenvalues 


+/- 15.384 
rad /sec 

+/- 2.3 1 lz 
rad / sec 

+/- 14.841 
rad/sec 

+/- 3.8141 
rad /sec 

+/- 1.930 
rad / sec 

States 

Mode 1 
eigenvectors 

Mode 2 
eigenvectors 

Mode 3 
eigenvectors 

Mode 4 
eigenvectors 

Mode 5 
eigenvectors 

^ y 

+/+ 1.000 

+/+ 1.000 

0 

0 

0 

n z 

0 

0 

+/+ 1.000 

0 

0 

9y 

+/- 0.065 

-/ -1-0.4331 

0 

0 

0 

e z 

0 

0 

+/- 0.067 

0 

0 

Vx 

-/- 0.043 

+/+ 0.598 

0 

0 

0 

Vy 

0 

0 

0 

0 

+/+ 1.000 

v z 

0 

0 

0 

+/+ 1.000 

0 

X 

-/+ 0.003 

— /+ 0.2591 

0 

0 

0 

y 

0 

0 

0 

0 

+/- 0.518 

z 

0 

. 

0 

.i 

0 

-/+ 0.2621 

0 


27 




Figure 1. Initial alignment of coordinate system for Large-Gap Magnetic Suspension System. 
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Figure 2. Open-loop modes for five-coil configuration. 
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Figure 5. Variation of equilibrium suspension currents with yaw for seven-coil configuration 
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Figure 6. Variation of equilibrium suspension currents with yaw for eight-coil configuration. 


S 




